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We suggest a variant of the recently proposed experiment for the generation of a new kind of 
Schrodinger-cat states, using two coupled parametric down-converter nonlinear crystals [F. De Mar- 
tini, Phys. Rev. Lett. 81, 2842 (1998)]. We study the parametric oscillator case and find that 
an entangled Schrodinger-cat type state of two cavities, whose mirrors are placed along the output 
beams of the nonlinear crystals, can be realized under suitable conditions. 
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I. INTRODUCTION 

Schrodinger-cat states |3,0] are most important in the 
domain of fundamental quantum mechanics, since the 
study of their progressive decoherence 0,0] would pro- 
vide a better understanding of the transition from the 
quantum to the classical world ]p[. However, due to their 
extreme sensitivity to the decoherence caused by the in- 
teraction with the environment, such linear superposi- 
tions of macroscopically distinguishable states are diffi- 
cult to produce and to observe ]p|fl] . In the last few years, 
a major effort in this field has led to the experimental 
production and detection of mesoscopic superpositions 
of distinct states, both in the context of the single- mode 
microwave cavities [^ and of the dynamics of the cen- 
ter of mass motion of a trapped ion B. On the other 
hand, entanglement has been widely recognized as one 
of the essential and most puzzling features of quantum 
mechanics M, in that it allows the existence of quan- 
tum correlated states of two noninteracting subsystems: 
Entangled states play a crucial role in the so called Ein- 
stein, Podolsky and Rosen (EPR) paradox [|j, and are 
essential in the rapidly growing field of quantum infor- 
mation, as they allow the feasibility of quantum state 
teleportation |^, quantum cryptography |1^, and quan- 
tum computation |11| . 



In two recent papers |12|, one of us has proposed an 
original scheme for the generation of a new kind of am- 
plified Schrodinger-cat type states. It is based on the new 
concept of quantum injection into an optical parametric 
amplifier (OPA) operating in entangled configuration. 

As a relevant variant and a natural extension of the 
above scheme, in the present work we analyze the case 
of the quantum injection in an optical parametric oscil- 
lator (OPO) in which two optical cavities are added to 
the OPA scheme considered in Ref. ]0]: refer to Fig. 0. 
Since the presence of the cavities leads to a large enhance- 



ment of the nonlinear (NL) parametric interaction, the 
number of the photon couples which are expected to be 
generated, in practical conditions, by the OPO scheme 
is far larger than in the amplifier condition: In addition, 
the generation of parametrically coupled quasi-coherent 
fields represents in this context an appealing perspective. 

The Schrodinger-cat state that has been put forward 
in Ref. ]|l^ and is being analyzed, in a more detailed 
fashion, in the present paper, is a superposition of two 
macroscopic states which are distinguished by their po- 
larization. It can be considered as a sort of amplified 
version of the polarization-entangled states which have 
been widely used in the last few years for the demon- 
stration of the violation of Bell's inequality ]0;n4|, of 
teleportation |^ , and for the generation of Greenberger- 
Horne-Zeihnger (GHZ) states |l^. 

The present paper is organized as follows: In Sec. || we 
briefly describe the process of type-II parametric down 
conversion, with an emphasis on the kind of entangled 
states usually produced in these experiments, and on the 
state we want to generate. In Sec. [II we outline the 
experimental apparatus needed for the realization of our 
scheme. We devote Sec. [V to the presentation of the 
dynamical time evolution of the density matrix and of 
the Wigner function in our system, and Sec. ^ to the 
discussion of the stability conditions for our paranret- 
ric oscillator. In Sec. ^ we set the initial conditions 
for the two coupled nonlinear crystals and the two cavi- 
ties, whereas the way in whic h the cat s tate i s produced 
is discussed in detail in Sec. VII. Sec. VIII is devoted 



to the presentation of the three methods we propose 
for d etecting the Schrodinger-cat state: photodetection 
(Sec. VIII A ) , measure ment of the second-order quantum 
coherence (Sec. VIII B ), and Wigner function reconstruc- 
tion (Sec. VHIC). We finally summarize and discuss our 
results in Sec. IX. The appendix is devoted to the devel- 
opment of the small interaction time approximation. 



II. ENTANGLEMENT GENERATING 
PARAMETRIC DOWN CONVERSION 

Let us first describe the kind of states commonly gen- 
erated in the experiments aimed at the violation of the 
Bell's inequalities. In these experiments the NL crystal 
(typically beta-barium-borate: BBO) is cut for Type II 
phase matching where the two down-converted photons 
are emitted into two cones, one "ordinary" polarized (o), 
the other "extraordinary" polarized (e). When the an- 
gle between the pump direction and the nonlinear crystal 
optical axis is sufficiently large |13|, the two cones mutu- 
ally intersect along two lines, lying on opposite sides of 
the pump beam direction. These ones identify the out- 
put modes of the parametric down conversion: kj (j — I, 
2). Therefore the field belonging to the modes kj can be 
simultaneously e- and o-polarized. In typical conditions, 
the output state of the emitted photon couple may be 
expressed by l9|JT^,|lF 



1 



lip) = -y= (|ei, 02) + e'"^|oi, 62)) 



(2.1) 



Since we have, for each couple, four degrees of freedom 
involved, i.e. 2 states of orthogonal linear polarization 



e, o for each mode kj, we can rewrite state (2.1) in the 
more precise form 

IV) = 4s (Il)le|l)2o|0)2e|0)lo + e'* 1 1) lo | l)2e |0) le |0)2o) , 

v2 

(2.2) 

which will be used in the following. 

The "Schrodinger-cat state" we want to generate is a 
sort of amplification of this state, that is, it may be ex- 
pressed in the form 

1 



|V')--^(|Ole|V^)2o|0)2e|0)lo 

+ e»^|Olo|V^'^)2e|0)le|0)2o) , 



(2.3) 



where \ip^) is a state with a large number of photons in 
some sense, and the states |0) are to be interpreted here 
as squeezed vacuum states. This kind of state is differ- 
ent from the traditional Schrodinger cat states discussed 
in the quantum optics literature Mn, where one has a 
single mode of the electromagnetic field in a superposi- 
tion of two macroscopic states with different phases of 
the field. The state ( |2.3| ) is a nonlocal superposition in 
which a macroscopic optical field is "localized" simulta- 
neously either in the e- or in the o-polarized mode. In 
other words it is a state more similar to nonlocal field 
states such as 



^(|a)i|0)2 + |0)i|a)2) 



(2.4) 



We shall present here an experimental scheme for the 
generation of a state which is actually a mixed state, 
but non etheless, has the same structure of the state of 
Eq. (2.3), that is, can be represented by the density op- 
erator 

P = 2 (P(^)le,20 «) p(0)lo,2e + P(0)le,2o «) p(A^)lo,2e 



+p(INT)ie,2o®p(INT')lo,2e 
+P(INT)L,2„ ® P(INT')L^2e) , 



(2.5) 



where p{N) is a two-mode mixed state with a large num- 
ber of photons, p(0) is a two-mode mixed state with a 
small number of photons and (o(INT) and /9(INT') are 
the interference terms. 



III. THE EXPERIMENTAL SCHEME 

We shall consider an experimental arrangement. Fig. [l|, 
based on the one proposed in Ref. |12] and similar to 
that adopted in Refs. [Q to show the realization of in- 
ducing coherence, without induced emission. Two down- 
conversion NL crystals, are arranged in such a way that 
the two corresponding idlers beams are aligned along a 
common direction fc2- Moreover, both idler beams and 
the signal beam of one NL crystal (with wave- vector ks) 
are placed within couples of mirrors. This scheme can be 
thought of to realize the coupling of two nondegenerate 
OPOs. The signal beam of the other crystal, emitted 
along the direction fci triggers the photodetector Di. 

The directions fci, ^2 and k^ are selected to realize for 
both NL crystals the Type-II phase matching described 
before. These beams are then associated with six modes, 
with annihilation operator aio, aie, 020, a2e, 030 and a^e- 
Note that the first two annihilation operators refer to 
traveling- waves, while the last four refer to cavity modes. 

The dynamics of the system is determined by the non- 
linear parametric interaction at each crystal and by the 
damping terms associated with losses and dissipation in- 
side the cavities p9] , as we shall see in the next section. 



IV. TIME EVOLUTION FOR THE DENSITY 
MATRIX AND THE WIGNER FUNCTION 



The partial Hamiltonian operators describin g th e uni- 
tary dynamics inside the crystals are given by |20 



where the field can be simultaneously in one cavity or in 
another cavity and whose generation is discussed in pTl . 



Hnl2 = i^X2(aL^L ^ 0,2oa3e) 
+ihx2{al^al^ - a2ea3o) 



(4.1a) 
(4.1b) 



where xi = £iX^^\ X2 = e2X^'^\ X^^^ is the second-order 
nonlinear susceptibiUty of the crystals, and e^ (i = 1, 2) 
is the pump intensity in crystals 1 and 2, respectively, 
which is assumed to be "classical" . 

Due to the explicit presence of dissipation in this prob- 
lem, one has to write the master equation for the reduced 
density matrix of the co mbine d syste m wh ich arises from 
the Hamiltonian terms ( 4.1a ) and ( 4.1b ) and from the 
damping terms 



A/3 = Hi{2aipa, 



?A 



a\aip 



.f.U. 



pa] a,) , 



(4.2) 



for i — 2e, 2o, 3e, 3o. Since the damping constants Ki are 
essentially connected to the transmittivity of the mir- 
rors, it is quite natural to assume K2e = 1^20 — 1^2 and 

Upon writing the full master equation for the total den- 
sity matrix p^ of the (six-mode) system, it appears clear 
that the dynamics of the six modes actually decouples 
into two independent dynamics for two groups of three 
modes. In fact, one has 



PT — ^le-2o-3ePT + ^lo~2e~3oPT 



(4.3) 



where 



real variables (or three complex variables). Its time evo- 
lution, upon evaluating the commutator and the damping 
terms and after some lengthy algebra, is described by the 
six-dimensional Fokker-Planck equation 

W{z,t) = 7,,— iz,Wiz,t)) + D,,^^-^Wiz,t) , 



dt 



dz. 



dzidzj 



(4.8) 



where the vector z — {xi,yi,X2,y2,X3,y3), the matrix 
D = diag(0, 0, K2/4, K2/4, K3/4, K3/4), and 



/ -XI \ 

xi 

-Xl K2 -X2 

Xl K2 X2 

-X2 K3 

V X2 Kg/ 



(4.9) 



The solution to Eq. ( [4.8|) can be written [g3[ as the 
integral 



W{z,t)^ d*z'W{z',0)T{z,z',t) 



(4.10) 



where 



^le-2o-3ePT — — ^[Hie-2o~3e, Pt] 



(4.4) T{z,z\t)^ 



+ K2(2d2oPTa2o - a2oa2oPT - PTa2o^2o) 
-I- K3(2a3ePTaL ~ 4e«3ePT " PTflL^Se) 



and 



Hle-2o-3e = «%l(aLa2o ~ "■lea2o) 

+i?iX2(aLaL ~a,2oa3e) 



(4.5) 



-^10-26-30 is identical to £ie-2o-3e up to the substitution 
e — > o and o —^ e. As a consequence, the complete time 
evolution will be of the form 



(27r)3 ^Deta(t) 



X exp 
G(t)= exp(-7i) 



1 



(4.11a) 



--{z-G{t)z'\a-\t)\z-G{t)z') 



and 



G* being the transposed of the matrix G. 



(4.11b) 



Z 

<j{t) ^ 2 I dTG{T)DG\T) , (4.11c) 



PT(i) 



=*g£i 



*Pt(0) 



(4.6) 



From Eq. ( |4.6| ) it is clear that if the initial condition is 
factorizcd, namely, if 



Pt(0) = pie-2o-3e(0) ® /Olo-2e-3o (0) 



(4.7) 



the state will remain factorized at all times, unless specif- 
ically designed conditional measurements pl[ are per- 
formed on the system (for example on the mode fci). 

Due to the decoupling between the le — 2o— 3e and the 
lo — 2e — 3o modes, we can simply restrict ourselves to 
the investigation of th e th ree-mo de p roblem described by 
the master equation (4.3) with (4.4), and we shall drop 
the subscript e and o when not needed. 

The Wigner function [|2) M^(a;i, yi, 0:2,^2,2:3, ys) = 
W{ai,a2,cx.2), with Ui = Xi + iyi {i = 1,2,3), resulting 
from this density matrix p will then be a function of six 



V. STABILITY 

The stability properties of the system are intimately 
connected to the threshold of the overall OPO consisting 
of NLl and NL2. Below threshold, the system is sta- 
ble and reaches a stationary state, since all eigenvalues 
of 7 have positive real parts. On the other hand, above 
threshold the system is unstable and its energy expo- 
nentially increases, because some eigenvalues of 7 have 
negative real part. 

This result can be easily checked in the case in which 
the parametric oscillator associated with NL2 is de- 
coupled from NLl (xi — 0): in this case modes k2 
and k^ decouple from mode fci, and we end up with a 
four-dimensional problem for the modes fc2 and k^, de- 
scribed by a Fokker-Planck equation of the same type as 
Eq. (Esl), but with 



«^2 





-X2 








K2 





X2 


-X2 





K3 








X2 





K3 



(5.1) 



and D = diag(K2/4, K2/4, K3/4, K3/4). In this case the 
(doubly degenerate) eigenvalues of 7 are 



K2 + K3 



± 



K2 - K3 



xi 



and the stability condition becomes 
X2 < '^sKs , 



(5.2) 



(5.3) 



which coincides with the customary threshold for the 
parametric oscillator pG]. However, if we turn on the 
first parametric amplifier {xi 7^ 0) then the problem 
turns from 4-dimensional to 6-dimensional, as we have 
seen: the eigenvalues of 7 change and it is in principle 
possible to change the threshold, i.e., the stability con- 
dition. As soon as xi 7^ 0; namely the first parametric 
amplifier is present, the system becomes unstable, inde- 
pendently on the values of X2^ 1^2, and K3. In fact, the 
eigenvalue equation for 7 is 

(A^ - (k2 + K3)A2 + (K2K3 -xl- X2)A + Ksxlf = . 

(5.4) 

As a consequence, we have three doubly degenerate eigen- 
values (Ai, A2, and A3). Since A1A2A3 — — K3X1, at least 
one of the Ai has a negative real part. 



VI. CHOICE OF THE INITIAL CONDITION 

We assume that at the beginning the first crystal is 
switched off (the pump strength ei =0). On the other 
hand, the second pump is on (e2 7^ 0) and the second 
parametric oscillator is in its equilibrium state below 
threshold. We therefore have a factorized initial state 



Pt(0) = Ple-2o-3e(0) «> Plo-2e-3o(0) 



(6.1) 



where 



Ple-2o-3e(0) = pio-2e-3o(0) = pi_2~3(0) 

= |0)ii(0|<^P2-3(0), (6.2) 

and ^2-3(0) is the equilibrium state of the oscillator be- 
low threshold. This can be easily determined upon con- 
sidering the limits 



lim G{t) = , lim cr(t) = cr(oo) , 

t — *CKj t — ^00 

and results in the following expression 
W{z,t = 00)^ I d'^z'W{z ',0)T{z,z', 00) 



(6.3) 



(6.4) 



1 



(27r)2 ^dctCT((X)) 



exp 



za ^(oo)z 



The equilibrium state is thus a Gaussian state in which 
the modes k2 and fca are correlated. 

The initial state ^2-3(0) is then given by the density 
matrix corresponding to the Wigncr function 



<t (0) 



2\V^_X| 

fc2 



X exp {-2 (x^ +yi + xl + yl - 2^x2X3 
+2^2/22/3)} , (6.5) 

where it is straightforward to realize that the modes 2 
and 3 are correlated. Moreover, it is not a pure state, 
because 

Tr(p2) = TT J dx2dy2dx3dy3 [W^tt'^W]' 



-'1-fl^Ki 



(6.6) 



as expected. 

The reduced density matrices of each mode are identi- 
cal and coincide with the thermal state 



xi 



^t'^(o) = - 1-f )-pr2|ani-f 



with an initial mean number of photons given by 



N 



xl 



2(fc^"Xi) ' 



(6.7) 



(6.8) 



which means that when the oscillator is initially suffi- 
ciently close to threshold the initial mean number of pho- 
tons in modes 2 and 3 within the cavities can be large. 



VII. GENERATION OF THE CAT STATE 



At time t — Q the first pump is turned on (ei 7^ 0): Also 
the first crystal begins to operate and the two groups of 
three modes start their joint evolution, according to 



PtH) 



-"=Vle-2o-3e(0) 
^=-^°Vlo-2e-3o(0) , 



(7.1) 



where the two factorized evolutions are identical because 
both the operator C and the initial condition are identi- 
cal in the two cases. As a consequence, we end up with 
two identical six-dimensional problems. 

The solution of Eq. ( [7.l| ) can be found as in Sec. f^by 
using the Wigner functions 



Wu3{z,t) = / d''z'Wi23{z',0)T{z,z',t) , (7.2) 



where the initial Wigner function Wi23( z^ ) correspond- 
ing to the initial density matrix [Eq. (|6.2[)] is given by 



Wi23{z,0) 



-) fl- 








exp < — 2 


2^2 + ^2 + 2:3 + yl 




-2 -r (3^32:2 + y3y2) 


} 


(7.3a) 


- ) Vd^ 


Ce^p{-2{z\C\z)} , 


(7.3b) 



with 



C 



and 






k 



1 





detC^(l-||' 



/I 














1 














1 














1 









_X2 

k 






X2 

k 






X2 

k 


1/ 



(7.4) 



(7.5) 



From Eq. (7.S) one can immediately recognize that 
since the initial state Wi23{z, 0) is Gaussian and the prop- 
agator T{z,z',t) is also Gaussian, the Wigner function 
Wi23{z, t) of the evolved state must remain Gaussian at 
all times. 

Upon integrating over d^z', Eq. {TJl) can be rewritten 
as 



W{z,t) ^ ^/i2*Mexp{-(zlB(i)|z)} , (7.6) 



where 



B{t) 



2a{t) 



G{t)C-^^G\t) 



(7.7) 



and G{t) and a{t ) are the six-dimensional matrices de- 
fined in Eqs. (Q, ( |4.11bD , and ( ^.llc| ). 

This Gaussian evolution holds for a short time only. 
As a matter of fact, one should distinguish between 
the mode along direction 1 and those along directions 
2 and 3: a!^ and Og denote creation of a photon in 
the stationary-wave modes within the cavities, whereas 
a[ denotes the creation of a photon in the traveling- 
wave mode along direction ki. Therefore the interaction 

J^NLi = i^Xi {O'l'A — 0,10-2) exists only for the time period 
during which this traveling wave mode 1 moves within the 
nonlinear crystal. In order to prepare the desired state for 
the modes 2 and 3, simultaneously taking full advantage 



of the degree of freedom represented by the traveling- 
wave mode 1, we perform a conditional |2l| ] measurement 
on direction 1, thereby conditioning the state of the four 
modes along directions 2 and 3 upon the detection of a 
photon along direction 1 polarized at 7r/4 with respect 
to the two output polarizations e and o, which are or- 
thogonal to each other. In this way we also post-select 
(along direction 2) the input state of the second crystal. 
The projection operator associated to such a conditional 
measurement is therefore given by 



PlL 



1 



{|l)lo|0>le + |0)i„|l)ie} 



X{lo(l|le(0|+lo(0|ie(l|} 



(7.8) 



As a consequence of this measurement (whose success 
probability amounts to 0.5) the state along direction 1 
and directions 2 and 3 factorizes: The state along direc- 
tion 1 is given by 



V2 



{|l)lo|0)le + |0>io|l)ie} , 



(7.9) 



which represents a photon polarized at 7r/4, whilst the 
conditional state for directions 2 and 3 is represented by 
the density matrix 

P2o-3e-2e-3o{t) « {lo(l|le(0| + lo (0| le (1 1 } Plo-2e-3o (*) 
(8)pie^2o-3e(i){|l)lo|0)le 

+ |0)io|l)ie} , (7.10) 



which can be rewritten as 

P2o-3e-2e-3o\^) ^ 



(1) (0) (0) (1) 

P2e-3oP2o-3e + P2e~~3oP2o-3e 



(int) (int) t (int) f (int) ,„ ... 

' /^2e-3oF2o-3e + /^2e-3oF2o-3e ; \'-^^) 



where 



P2-3 = l(l|Pl-2-3(t)|l)l, 
4°^3 = l(0|Pl-2-3(t)|0)l, 
= l(l|pi-2-3(i)|0)l. 



(int) 
P2-3 



(7.12a) 
(7.12b) 
(7.12c) 



The state o f Eq . ( [7.11 ) is of the same form of the desired 
state, Eq. ( |2.5D , and is a linear superposition of distin- 
guishable states, as long as P2-3 is well distinguishable 

from p-^lj^. 

It should als o be emphasized at this stage that the 
density matrix ( [7.11 ) directly corresponds to the Wigner 
function, Eq. (2) of Ref. |1^, obtained in the OPA case. 
The similarity between the OPO and the OPA configu- 
rations is better brought about in the limit of small in- 
teraction times (see the appendix, where it is also shown 
that — in this limit — many of our results are very simi- 
lar to those obtained in the OPA case ||l^). Roughly 
speaking, one should recover the OPA results from the 
OPO ones in the limit k -^ 00, since this condition means 



absence of cavity mirrors. However, this correspondence 
does not hold exactly because the initial state in the OPO 
case (the state present in the cavity at t = 0, when the 
first nonlinear crystal is switched on) is slightly different. 
This fact explains the differences between the OPA and 
the OPO, which result in a far larger effective number of 
photons in the latter case. 



VIII. DETECTION OF THE CAT STATE 

How can we probe the quantum state produced in this 
parametric-oscillator entangled configuration, and prove 
that it actually represents a Schrodinger-cat state? In or- 
der to do this, one has to independently show that i) the 
state is indeed made out of two macroscopically distinct 
components, that ii) these two components exhibit quan- 
tum interference, so that the state can be considered as a 
true linear superposition rather than a statistical mixture, 
and that iii) the "separation" between the two compo- 
nents scales with a macroscopic or mesoscopic parameter, 
usually the number of photons. To achieve this goal, we 
propose three different and independent methods — which 
can be used either alternatively or simultaneously — as we 
shall explain in detail in the next three subsections. 



A. Photodetection 



P2e 



-Lr2-3P2-3 -Lr2-3P2-3. 



(8.4) 



with an identical form for the reduced state p2o- The 
reduced density matrices p^e and p^o can be determined 
in a similar way. 

From Eq. {^A) it is immediate to recognize that the 
reduced state of the mode 2e is given by the sum of two 
density matrices, conditioned upon the detection of one 
photon and of zero photons in the mode la (or, more pre- 
cisely, one photon in the mode le), respectively. There- 
fore the two terms of the reduced density matrix can be 
experimentally obtained by rotating the polarizer in front 
of the detector Di located along the direction ki: When 
the polarizer is vertical (mode le), we have zero photons 
in the mode lo, and only the second term of the sum in 
the right-hand side (rhs) of Eq. (B^) is realized. On the 
contrary, if the polarizer is set horizontally (mode lo), 
one detects one photon in the mode lo, projecting the 
resulting density mat rix for the mode 2e onto the second 
term in the sum (B.4). However, both terms are present 



when the polarizer is set at 45°. An experimentalist could 
then take advantage of this property to test the presence 
of the two component states: the distinction between the 
two states in the superposition can be made via photon 
number measurements, yielding the probability distribu- 
tion P(n2e). In fact, one has 



Let us employ photon number measurements for the 
modes along direction 2, thereby collecting the photon- 
number distributions P{n2o) and P{n2e)- We therefore 
consider the reduced density matr ix ob tained by perform- 
ing the trace on the state of Eq. ( |7.1l| ), that is. 



P2e — Tr2o-3e-3o {P2e-3o-2o-3e(^)} 

= I {Tr3„ 



P2e-3o 



Trso 



P2e-3o 



Tr23 
Tr23 



(0) 

P2-3 



pi% 



}[-(I) 



where P(t: /A) is the probability of finding one photon 
with polarization at 7r/4, that is. 



P^Plo~2e-?,o{t) 



^ ( J ) — Trio_ie_2e-3o-2o-3e 



® Ple-2o-3e{t)] 



Tr9_.c| 



f^3 



Tr2-3 



(0) 

P2-3 



represents the probability of the conditional measure- 
ment generatin g the desired cat state. The interference 
term s in Eq. ( [7.11 ) obviously give no contribution to 
Eq. (8T), since 



Tr^ 



2-3 



(int) 
P2-3 



= Tr 



2-3 [1 



(l|pi-2-3(t)|0)i]=0. (8.3) 



Combining Eqs. (pTTl) and (8.2), one obtains for the re- 
duced state 



1 



P{n2e) ^ 7^{PYi{n2e) + PY{n2e)) 



(8.5) 



where PH("-2e) [^v("-2e)] is the probability distribution 
obtained when the polarized is set horizontally (verti- 
cally). The results an experimentalist would obtain with 
a simple photodetection in these two situations are shown 
in Fig. I (a) and (b) , together with the probability distri- 
bution ( S.5I) one would obtain when the polarizer is set 
at 45° [Fig. 1(c)] . 

In this way we have verified the existence of two dis- 
tinct components in the state (3.4). But how can we 
be sure that these two components form a quantum su- 
perposition and not just a classical mixture? To answer 
this question, one has to perform a measurement able to 
distinguish the "cat state" 



P2o-3e-2e-3o(^) 



2Tr2_3P2^-3Tr2-3P2"-3 



(1) (0) , (0) (1) 

P2e-3oP2o-3e + P2e-3oP2o-3e 

(int) (int) f , (int) f (int) 
P2e-3oP2o-3e + P2e-3oP2o-3c 



(8.6) 



from the corresponding statistical mixture 



P2o-3e-2e-3o(^) ~ 



2Tr2-3P2^-3Tr2-3P2-3 



(1) (0) I (0) (1) 

^ P2e-3oP2o-3e + P2e-3oP2o-3e ' 



(8.7) 



which does not exhibit any interference. 

In order to reach this goal, we perform an interfer- 
ence experiment, involving the modes along direction ^2 
only, using a detection system similar to the one pro- 
posed in Ref. [|l^, as schematically described in Fig. |^. 
The measured quantity is given by the photocounts at 
the detector Dc, as a function of the variable phase (p. 
The annihilation operator c corresponding to the mode 
traveling to the detector Dc can be written in terms of 
the annihilation operators of the modes 2e and 2o as 

c = -^ (a2o + e"f'a2e) , (8.8) 

so that the operator number of photons for the mode c 
will be given by 



c'c - 



W t t 

2 [a-2oO'2o + fla 



2ea2e + c"^ al^a2e + & ^"^ a\^a2o 



(8.9) 



In order to be able to distinguish between the superposi- 
tion state and the mixture, the expectation value 



(ctc)cat = Tr [c'^cplt 
has to be different from 



(ctc)„,i, = Tr [ctcp™!! 



=-3o(i)] 



3e-2e-3i 



oit)] 



(8.10) 



.11) 



It is then clear that this interference experiment can an- 
swer our question whenever the contributions of the off- 
diagonal terms Tr[c^cp2r-3oP2"-3e] ^'^'^ i^^ complex con- 

jugate T^nc^ cp\e-3oP2o-3,e\ ^re nonzero. 

Let us start by evaluat uig t he contribution of the diag- 
onal terms, namely, Eq. ( B.ll ). After explicit integration 
of the corresponding Wigncr function, it is easy to prove 
that the phase-dependent terms [the third and the fourth 
term in Eq. ( ^.9| )] vanis h when one computes the expec- 
tation value, Eq. (8.11). Therefore the diagonal terms 



yield a phase (0)-independent contribution given by 

(c'''c)niix = ^ ( (a2o'^2o)mix + (a2e"2e)mix 



(n2o><^) + (n2o) 



(0) 



(n2e>(^'+(n2e) 



(0) 



("2 



,(1) 



{n2 



,(0) 



where 



Tr2-3 



^2 



' (i) t ' 
P2-3«2«2 



Tr2-3 P2-3 



.12a) 
.12b) 

8.13) 



(« — 0, 1) is the mean photon num ber in one of the 
two diagonal states in Eqs. (S.6) and (8.7). In the small 



interaction-time limit, which is very well justified in the 
present case (see appendix), 1 3> fci, xit, X2^i we have 



p(%^P2-3(0), 

P2-3 o^ 4^2-3(0)02 



.14a) 
.14b) 



As a consequence, the two expectation values in the rhs 
of Eq. ( ^.121: ) can be explicitly evaluated and are given 

by 



\(o) 



(jl2)'"' = N 

(n2)(i) =2iV+l, 



^'S 



.15a) 
.15b) 



where N is the initial mean photon number in the cavity. 
In conclusion, the diagonal contribution to the expecta- 
tion value in Eq. (S.IOD amounts to 



{c'c 



l + 3iV 



(8.16) 



which is indeed ^-independent as expected. 

We t urn now our attention to the off- diagonal terms 
in Eq. ( |8.6|) , which are absent in Eq. (8/7). First we note 
that the expectation values of the number operators rela- 
tive to the two polarizations in mode 2 computed on the 
off-diagonal terms vanish, i.e., 

(«La2o)o-d = (aL«2o)o-d = . (8.17) 

On the other hand, the third and the fourth term in the 



rhs of Eq. (8^) give to the expectation value on the off- 
diagonal terms the contributions 



s"^(4oa2e)o-d = |2Tr2-3 



P2-3 



Tro- 



2-3 



F2-3 



>a2c 



\(int)/ t \(int)t 



{a2t 



^20/ 



v(int)t/ 



^20/ 



v(int) 



-icp 



(4e«2o)o-d = (e*"^(4oa2e)o-dj , 



where 



("2) 



(int) 



Tr^ 



2-3 



(int) 
^2-3^2 



(8.18a) 
(8.18b) 

(8.19) 



These contributions are generally different from zero, and 
this observation is sufficient to reach the conclusion that 
the proposed interference experiment is able to distin- 
guish the cat state from the corresponding mixture. 

We are able to evaluate these off-diagonal terms in the 
small interaction-time limit developed in the appendix: 
At the lowest order in xit, X2t, and kt, we have 



(™t) t l'r\\ 

P2-3 -02^2-3(0) 

and therefore, using Eq. ( ^.IQD , 



.20) 



(a2e)('"')=Xii(A^ + l) , 


(8.21a) 


(a2e)('"*)^=Xii(«L>=0, 


(8.21b) 


(4o)^'"'^-xi^(4o> = o, 


(8.21c) 


(4o>^'"*^^ = xi^(^ + i) . 


(8.21d) 


On the other hand, 




P2'-3 - X?t'4P2-3(0)a2 , 


(8.22a) 


p("J3^P2-3(0), 


(8.22b) 


which yield, respectively. 




Tr2-3P^'23 = Xlt^ (N + l) , 


(8.23a) 


Tr2-3P2"-3 - 1 ' 


(8.23b) 



and, finally, 



J<P 



(4oa2e)o-d = 



N + l 



,i<P 



(8.24) 



In conclusion, considering the off-diagonal contribution, 
Eq. (|8.10) can be rewritten as 



(c'l'c) cat 



iC'C 



N + l 



(8.25) 



It is then clear that the photocounts at the detector Dc 
exhibit interference fringes as a function of the variable 
phase (j), if and only if the state (8^) is a true linear su- 
perposition and not just a statistical mixture of the two 
macroscopic components. The visibility of such interfer- 
ence fringes is given by 



V = 



1 + N 



(8.26) 



l + 3iV 
and has therefore the lower bound 1/3 for N ^ oo. 

B. Correlation functions 

Our aim in this subsection is to compute the first- and 
second-order correlation functions relative to our output 
modes, in order to make an independent test of the pres- 
ence of quantum coherence in our system. We keep in 
mind Eq] that a manifestation of quantum coherence at 
second order is subpoissonian statistics, i.e., 



G(2)(o)< [g(i)(0) 



.27) 



where G^^''(0) and G'^^\0) are, respectively, the first- and 
second-order correlation functions. 

Let us consider the same experimental apparatus we 
have proposed for the detection of interference (see 
Fig. y). We take now into account both output ports 
c and d of the polarizing beam splitter, with annihilation 
operators 



71^""'' 



i<t> \ 

e ^a2e) , 



d = — = (a2o - e*"*a2e) , 



;.28a) 
.28b) 



and evaluate the correlation functions (c^cc^c), (d'^dd^'d), 
and (c^cd^d), where c^c is given by Eq. (|8.9|), and 



d^d = - (al^a2o + al^a2e - e'"^4oa2e - e *"*4e"2o 

(8.29) 

We shall evaluate the functions (c'c)'^, (d^d)'^, and 
c^cd'^d in the small-time approximation limit (see ap- 
pendix), in which 



P2e-3o~2o~3e{t) OC (o^^ + a2oU P2e~3o{0)p2o-3e{0) 

X (a2e + a2o) , (8.30) 

where P2o-3e(0) is the Gaussian state described by the 
Wigner function W^^^{0) of Eq. {51), for which the 
Wigner function corresponding to the reduced density 
matrix of mode 2 alone is given by Eq. (6.7), that rep- 
resents a thermal s tate w ith a mean number of photons 
given by N of Eq. (B.15a). 

Upon evaluating all the required expectation values, 
we obtain 



{{c'c)') = 



{{duy 



16N^ + UN + 2 -h 2 cos {AN^ + 5N + 1) 
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(8.31a) 

167V2 + UN + 2 - 2 cos (4N^ + 5N+1) 



{c^cdU) ^N{2N +1) 



.31b) 
!.31c) 



From Eqs. ( 8.3l| ) it is clear that the visibility of the 
fringes in {{acy) and {{d^d)^) is given by 



V ^ 



4N^ 



57V +1 



8iV^ 



77V +1 



.32) 



and monotonically decreases from V = 1 (for A^ = 0) to 
F = l/2 (for TV ^(X)). 

Finally, considering the field at the output port c, the 
first- and second-order correlation functions for mode 2 
can be written as 



G«(0) = (ctc) 



l + 3N + {l + N)cos(l) 



G^^\0) = (c^c^cc) = {{c^c)^) ~ (c^c) 
^2N[l + 2N+ (iV+l)cos< 



(8.33a) 



.33b) 



respectively. It should be noted that these results map 
into the corresponding ones obtained in Ref. ^^ for the 



OPA case upon a redefinition of the phase angles. By 
comparing [G^^^(O)]^ and G^^^(O) it is possible to see that 
G^^\0) < [G(^^(0)]^ only at low mean photon number, as 
it could have been easily expected. The best situation is 
obtained when (h — 0, in which case 



G(2)(0) 



(1 + 2iV) 
2N (3iV - 



(8.34a) 
(8.34b) 



and the condition for quantum coherence at second order 
is reached when N < l/\/2. On the other hand, when 

= TT, [G(i)(0)]^ = ^^ G(2)(0) = 2iV2, and therefore 
G(2)(0) is always larger than [G^^\0)]^. 



C. Wigner function 

The aim of the present section is to provide a means 
to represen t the essential features of the Schrodinger-cat 
which "lives" 



state, Eq. (7.11), which "lives" in a 8-dimensional phase 
space, in the more customary 2-dimcnsional phase space, 
in order to make a comparison with the more conven- 
tional cat states [g|,|| . Let us start from Eq. ( 7.11 ) which 
we rewrite here for convenience 



P2o-3e-2e-3o(^) ^^ 



(1) (0) , (0) (1) 

P2e-3oP2o-3e + P2e-3oP2o-3e 

(int) (int)t (int)t (int) 1 /o oc\ 

P2e-3oP2o-3e + P2e-3o^2o-3e ■ [O-OO) 



The Wigner function representation of the density ma- 
trix ( p.35| ) would of course reflect its characteristic 
Schrodinger-cat properties. However, in order to better 
understand the nature of this state, it would be interest- 
ing and desirable to see whether it is possible to find dif- 
ferent optical modes in whose terms the state (and there- 
fore the Wigner function) may be rewritten in a simpler 
form. Our key idea is then to look for linear combinations 
of mode operators (which can easily be realized with lin- 
ear elements: polarizers and beam-splitters) such as to 
factorize the state ( p.35| ) in smaller subspaces. 

We first perform a transformation which changes the 
horizontally and vertically polarized modes into the 45°- 
polarized ones, namely, 



0-2e + 0-20 a,2e — 0'2. 

a+45,2 — j^ — , a-45,2 



x/2 

a+45,3 — y= — , a_45,3 



V2 

0-Ze — 0-3.0 



V2 



;.36a) 
.36b) 



and the corresponding expressions for mode fci and for 
the creation operators. In t erms of the se new operators, 
ffNLi and iJNL2 [Eqs. ( 4.1a ) and ( 4.1b| )] can be rewritten 
as 



-^NLl = iTT'Xl (^^+45, 2 '^-1^45 4 ~ a+45, 20+45,1 

-ihxi [0-45,20-454 ~ a-45,ia-45,2J , (8.37a) 

^NL2 = if>'X2 (0+45,2^+45.3 ~ «+45,2a+45,3 ) 

-1^X2 (aL45 20-45.3 ^"-45,20-45,3) ■ (8.37b) 



We have already assumed [Sec. |V| that the cavity decay 
rates k^ do not depend on the polarization. This in turn 
means that k+45_2 = ^^-45, 2 = 1^2 and K+45,3 = K-45,3 = 
K3, and therefore we have that for the ±45°-polarized 
modes we have the same evolution equation as that for 
the original modes (except for a minus sign). Conse- 
quently, it is possible to repeat all the same arguments 
as before [Sees. IV and VI] |. In particular, the modes 
a+45^1, a+45_2, and 0+45^3 are decoupled from their or- 
thogonal counterparts 0-45,1, a-45,2, and 0-45,3, and the 
evolution equation may be rewritten as 



PT{t) 



e^+«*e^- 



P+45,l:+45,2;+45,3;-45,l:-45,2:-45,3(0) . 

(8.38) 



38) the initial condition is given in the same 



In Eq. 

way by 



P+45,l;+45,2;+45,3(0) — |0) 



+45,1 \ 



P+45,2;+45,3(0) I 



(8.39) 



where p+45 2+45 3(0) is the equilibrium state below 
threshold of the parametric oscillator when NLl is turned 
off, and the same initial condition holds for the —45°- 
polarized modes. As a consequence, the same Gaus- 
sian evolution we have found in Sec. IV holds. The 



only difference is that now the conditional measurement 
is simply a projection onto the state |l)+45,i, i.e., the 
one-photon state for the 0+45,1 mode, while the —45°- 
polarized modes remain decoupled from the orthogonal 
ones. 

The cat state after the conditional detection of the 
n = 1 photon for the +45, 1 mode is then written in the 
following way 



(1) 



p^ ex -45,l(0|+45,l(l|P+45,l;+45,2;+45,3(*)P-45,l;-45,2;-45,3|l)+45,l|0)-45,l = P+45 



2:+45,3 



(0) 
' P-45,2:-45,3 



(8.40) 



where p^_c. and p\l^ a re again given by the expres- 



sions (7.12a) and ( 7.12b| ). It should be noted that, using 
these new ±45°-polarized modes, one gets a complete 
factorization of the — 45°-polarized modes, which are not 
affected by the quantum injection process induced by the 
conditional measurement. The — 45°-polarized modes are 
not "interesting" , in the sense that all the "cat" proper- 
ties of the state ( 8.40| ) are contained in p^45 2+45 31 ^-iid 
therefore we shall neglect them from now on. We are 
then left with the state pVis 2-+45 3) which is an entan- 
gled state of the modes -1-45, 2 and -1-45, 3. 

As the second step of our procedure aimed at the fur- 
ther simplification of the original 8-dimensional Wigner 
function, we consider the transformation 



which is suggested by the int eractio n term in Eq. 
In terms of d+ and (i_, Eq. ( 8.37b| ) becomes 



.37b) 



H, 



NL2 



ifl 



^(^n 



^dl' _ 4) _ ,n^ (dtj _ ^2_^ 



..42) 



and the two modes (i+ and d- are squeezed by the nonlin- 
ear crystal. These modes can be experimentally realized 
outside the cavity for example with two PBS and a 50%- 
50% BS, as schematically described in Fig. ^. The state 
of these two modes can be represented by the Wigner 
function 



Wi 



f Xd, + Xd_ Vd+ + 2/d_ Xd. - Xd_ Vd+ - yd- 



d+ = 



a+45,2 + a+45,3 



V2 



d- = 



a+45,2 — a-l-45,3 



V2 



,.41) 



V V2 ' a/2 ' V2 
where [see Eqs. {^) and ( [TT] )] 



V2 



Wi{x2,y2,X3,y3) ex / dxidyiWn=iixi,yi) x 



VdetB 



-{z\B\z) 



(8.43) 



(8.44) 



What is the nature of this state? In order to answer this question, we are naturally guided by two different 
approaches: i) the study of the OPA case [|l2[ and ii) the use of the small-time limit xii, X2t, nt <^ 1 we have already 
considered in Sec. VIII A and worked out in the appendix. In the OPA case [Oj the output state at time t is given by 



IV'(t)) = ^e^^*("2e4o-'^2ea3o) + X2t(aLaL-''2oa3e) ^ L^^ + flH |0) 



which can be rewritten in terms of the ±45°-polarized modes as 



X2t(al45,2°-45,3-<*-«.2a-45,3)|r)\ g^ X2t(Q^45_2oJ^45_3-a + 45,2a + 45,3) t \n\ — ^/.C) „/,(!) 

\^/'i^f^ T , T , "+45,2l^/ — ''^-45,2:-45,3'/^+45,2;+45,3 



Neglecting the factorized state V'-is 2-45 31 ^^"^ using the d± modes, we have 



V^+45,2;+45,3 



g22£(df-4)g-2^(dL^-d^_) 



'di+d'' 



V2 



|0) 



X2t ^ 

2 ' 



X2t 



f,0 



X2t 



(8.45) 
(8.46) 

(8.47a) 
(8.47b) 



which is an entangled superposition of the squeezed one- 
photon and vacuum states of the modes d+ and d- . It 
is quite clear now that if we want to "isolate" one mode, 
say, the ^4. mode, we need a second conditional measure- 
ment on the mode d_, e.g., a projection onto the state 



\v)d_ =a|0)rf„-f/3|l)rf_ 



(8.48) 



Such a conditional measurement could be performed, for 
example, by sending a two-level atom — resonant with the 
atomic transition — through the cavity, and eventually 
post-selecting its internal state in a corresponding super- 
position of its ground and excited states. The conditional 
state, provided the measurement has given a successful 
result, would then read as 



l«--«°*lf>«'* + S^,*'2 



^,0>i, , (8,49) 



We can reach a similar conclusion also by analyzing 
the OPO case using the very well justified small-time ap- 
proximation (see appendix) in the limit xit, X2t, Kt <ti 1, 
applied to the modes 4-45,1, -f45, 2, and -1-45,3. We 
have, at the lowest order in xit. 



,(1) 



P+45,2;+45,3 - +454(l|Pl-2-3|l 



2-31-1-/ +45,1 



^ a+45.2P2-3(0)a+45,2 , 



(8.50a) 
(8.50b) 



where the initial density matrix ^2-3 (0) is the state de- 
scribed by the Wigner function (|6.5[). If we now write 



10 



Eq. (6.5) in terms of the new variables corresponding to 
the modes d+ and (i_, namely, 



Xd+ 



X2 +XJ, 

V2 



Xd- 



X2 - X^ 



^/2 



.51a) 



Vd. 



we obtain 



y2 + y3 

V2 ' 



Xd- = 



?/ 2 -y3 

V2 



k: 



d_ 



<"<" 



^exp 



-2 1 



X2 



O-f) 



yi 



1 ^ exp 



-2 1 



^).L-2(1 



X2\ 2 



The initial states for the modes d± are generalized Gaussian states |^^ , of the kind 

p± (X exp ~nd]j.d± TO^cfj. m±d]j_ 

with 



log 



1 _ X2 
J^ ..9 



1 ^2 



1 X2 

^ ..9 , 



and 



, X2 

m± = ± — n 



Since the initial state factorizes, we have 

Pd+d- ^ a+45,2Pd+iO)pd_iO)a+45,2 oc {dl_ + dl)pd+iO)pd_iO){d+ + d-) 
= (4pd+ (0)d+) ® pd_ (0) + pd^ (0) ® {dlpd_ (O)dL) 

+ (4Pd+(0)) ® {pdAO)d-) + {Pd+{0)d+) ® (dLpd_(0)) , 



.51b) 



(8.52) 
(8.53) 

(8.54a) 

(8.54b) 

(8.55a) 
(8.55b) 



which is a mixed st ate a nalogous to the pure state ( B.47|) obtained in the OPA case. Its Wigner function can be 
calculated from Eq. (^.5E ) and is given by 



W {xd_i_,yd+,Xd_,yd-) = 



i^'^r 



'' 2^H 



-2 4 



xl 



K..L)(2-f)V(.^...g(2.^)^-2 

e-2(i-¥)K+i.L)-2(i+^)(=^L+v^J 



2 ; \Xd^Xd_ +yd+yd-) 



(8.56) 



Two important features should be noted within the form 
of this Wigner function: i) the interference term (the last 
term in the square brackets) decreases when the number 
of photons in the initial state increases. This behavior 
is governe d by the factor 4 — "x^/ k^ and by the fact that 
[see Eqs. (IS) and ( B.15 )] N ^ oo when X2/K -^ 1. ii) 
The Wigner function is negative around the origin and 
its negativity scales to zero as the initial mean photon 
number N ^ 00. In fact, 



1^(0,0,0,0) 



('-#)' 



2-# 



7r2 {2N + 1){N + 1) 



.57) 



We have already seen the same scaling behavior of quan- 
tum properties with A^ ^ 00 in the calculation of the 
second-order correlation function G*^^-* : This is one of the 
desired properties of a Schrodinger-cat state, as we have 
emphasized at the beginning of this section. 
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Again, Eq. (B.56) bears a remarkable similarity with 
the corresponding result obtained in Ref. [|2| for the OPA 
configuration, in the limits k ^ oo and of small interac- 
tion times. The main advantage of the OPO is given by 
the l arge r effective number of photons per mode N [see 
Eq. (Q] with respect to the sinh^ xt of the OPA ^. 

We have therefore learnt that in order to obtain a one- 
mode state which embodies all the relevant features of 
the original four-mode cat state one has to perform a 
conditional measurement on the mode d-. When this is 
successfully done, the final conditioned state of the mode 
d+ alone is described by the Wigner function 



M/(<5+) (XTT / d'S-W^\o)+p\i)i6-)Wi6+,S^) , (8.58) 
where 



W: 



(xd+ + Xd_ yd+ + yd_ Xd+ - xd_ yd+ - Vd_ 



'\^/f 



V2 



V2 



V2 



(8.59) 



is the Wign er function [see Eqs. ( ^.43[ ) and ( p.44[ )] of the 



state 



.55|), and 



W^a|0>+/3|1>('5- 



2 r 



1 + 4|/3|^(^|<5|^- 

-4Re((5)Re(a/5*) + 4Im((5)Im(a*/3) 

(8.60) 



is the Wigner function of the state onto which t he co n- 
ditional measurement projects the mode d- [Eq. ( |8.48 )] 
According to the small-time limit approximation (s ee ap - 
pendix) the explicit form of the W i gner function 



can be derived from Eqs. (^.55 )-( p.6C| ) and, after 
lengthy calculation, reads as 



5? 



X2 



W{xd^ , 2/<i+) (X exp -2 1 - ^ x^ - 2 1 + ^ yj 



'Vd~ 



X2 



- 1 



X2 




+ 2Re (a* (3 \xd+ (2 - — ) 



iVd. 



X2 



2-^ 



(8.61) 



which is in very good agreement with the numeri- 
cally computed exact one. As desired, the value of 
W{xdj^,yd+) at the origin may also be negative (depend- 
ing on the parameters a and (3 specifying the conditional 
measurement), reflecting the quantum pr oper ties of the 
original 4-dimensional Wigner function (S.56). Explic- 
itly, one has 



W{xd+ =0,2/d+ =0) 



X2 



.62) 



In 



These results are graphically shown in Figs. 0-u 
Figs. ^ and g we have plotted the Wigner function ( 8.61 ) 



for two values of A^. In both figures two different view- 
points have been selected for the tridimensional plots, 
in order to display most clearly the quantum superpo- 
sition character of our cat-like state. In particular, one 
should note that in both cases the Wigner function is 
negative around the origin. However, a comparison be- 
tween Fig. Q and Fig. ^ shows that even though the two 
Gaussian peaks are better separated for a larger number 
of photons, the negativity of the Wigner function tends 
to disappear as soon as the initial number of photons in- 
creases, as expected. This behavior is further confirmed 
by the inspection of the corre spond ing marginal distribu- 
tions of the Wigner function ( ^.61 ), shown in Figs, g and 
^: P{x) = P{xd^) displays a larger separation between 
the peaks as the initial mean photon number (rt) = N in- 



creases. On the other hand, P{y) = P{yd+) displays the 
interference between the two macroscopic components, 
which tends to be washed out when the number of pho- 
tons increases. In fact, for N = 14.94, the interference 
fringes have already disappeared. 



IX. DISCUSSION AND CONCLUSIONS 

In this paper we have considered the generation of 
entangled Schrodinger-cat states in an optical paramet- 
ric oscillator, as a relevant variant of the original pro- 
posal [|2| which instead had considered the amplifier 
case. In these works, the central point (both concep- 
tually and experimentally) is the quantum injection 112] 
of the second nonlinear crystal with the output of the 
first parametric medium. In the present paper, we have 
computed the time evolution for the electromagnetic-field 
and chosen the initial condition needed for the generation 
of the desired cat state. Such a state, however, lives in 
a eight-dimensional phase-space: therefore we have pro- 
posed three methods which are able to prove that it is 
an actual Schrodinger-cat state: direct photodetection, 
measurement of the correlation functions, and measure- 
ment of the Wigner function. Our calculations show that 
the state produced in this way has indeed two macro- 
scopic (mesoscopic) components — which are macroscopi- 
cally (mesoscopically) distinguishable — and that they are 
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in a coherent superposition (and not just in a statistical 
mixture), i.e. they display quantum interference. 

A comparison with the performance of the correspond- 
ing OPA scheme P2| is in order here. First, the OPO 
has a larger conversion efficiency due to the enhance- 
ment factor of the parametric interaction, given by the 
presence of the cavities. This leads to a larger number 
of photon couples with the same pump power. Second, 
our Schrodinger-cat state is confined in the cavities, con- 
trarily to what happens in the OPA case, where it is a 
traveling wave. However, the price one has to pay in 
order to have these advantages, is given by the unavoid- 
able cavity losses, that tend to destroy the coherence of 
the state when the number N of initial photons tends 
to infinity. Such a phenomenon — decoherence M-m — is 
visualized by the progressive disappearance of the inter- 
ference fringes and of the negativity of the Wigner func- 
tion when N increases. It is then clear that one has to 
consider a trade-off condition between the enhancement 
factor (a large N) and the losses (a low k). This may lead 
to a comparison between the performances of the OPO 
and the OPA ||lj]: in particular, our OPO configuration 
is preferable when the mean number of initial photons N 
[see Eq. (|6.8|) ] is larger than the corresponding parameter 



(sinh^ X* Til) of the OPA. 

In conclusion, we think that an experiment along the 
lines outlined in this paper and in |l2] — which is real- 
izable using presently available technology — is a promis- 
ing candidate for producing entangled superpositions of 
macroscopically distinct quantum states. 
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where Lk is the crystal length, n its refraction index, and 
c is the speed of light in vacuum. On the other hand, 
for an average pump power P ~ 300 mW, the coupling 
strength is of the order of X2 — 6 • 10® Hz. In order to ob- 
tain "macroscopic" states, one needs a quite large initial 
mean number of photons in the parametric oscillator be- 
low threshold. This fixes the damping rates K2 ^ Kg = k 
to be slightly larger than X2, since, from Eq. (6.8), we 
have 



xi 



1 



1 

2iV 



(A2) 



Therefore, we have k ~ 6 • 10® Hz, too. Since the wave- 
length of the photon is A ~ 7.3 • 10~^ cm, this amounts 
to having a standard cavity, with a quality factor 



2nc 
Xk 



W 



(A3) 



On the other hand, xi will be of the order of X2- In 
summary, we have 



X2t^Xit- i^i- 10" 



(A4) 



From Eqs. (|I^-|45|), (|t|), and (U), one has, for the 
time evolution of the combined density matrix. 



Pl23(0 



^^*V23(0)|0)i(0| 
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Since £123 oc k, xii X2j it is appropriate to expand the 
exponential e^^^^* in power series up to second order in 
i^t, Xit, X2t, yielding 
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1 + ^123* + ■;^-^123^ 



(A6) 
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APPENDIX: 

The fact that the time t during which we have the in- 
teraction within the first nonlinear crystal is very short is 
of fundamental importance, and it allows an immediate 
description of the experiment. To bring out this most 
clearly, we develop an approximate treatment, which is 
however justified by the actual experimental values re- 
ported in Ref. |l|. 

The interaction time t, which is the time of flight of the 
photon within the first nonlinear crystal NLl, is given by 



and 



t = 



Li,n 



10"" sec 



(Al) 



'C123P = C23P + X1 



alal — aia2,p 



(A7) 



where £23 is that part of the Liouvillian which only acts 
on the modes fc2 and fcs, as given by Eqs. (4^) and (|4.5|). 



It is possible in this way to determine the conditional 
states P23 ' P23 ' ''^^'^ the interference terms in Eqs. ( |7.1l| - 
7.12c|). We compute first 



Pi% 



4P2-3(0)a2 
Tr(4p2-3(0)a2) 



(A8) 



The properties of this state are usually characterized by 
measuring the photon-number distribution of the mode 
2 along direction 2. We have therefore to perform the 
trace over the mode 3 in Eq. (A8), obtaining 



(l)iod rj^ 

p\ = Tr 



^{p'^1.) 



Q|(Tr3p2^3(0))a2 
Tr(a^P2-3(0)a2) 



(A9) 



We already know that Tr3P2-3(0) is a thermal state with 
a m ean number of photons given by N [see Eqs. (p^), 
3), and (8.15a)], i.e., 



13 



Tr3P2-3(0) = ^fY 



N 



■N 



\n){n\ 



1 



1 + A^ 



[see Eq. (3/7)] and consequently 

oo 



n=0 



Pnin) = n 



N 



l + N 



1 



(1+iV) 



2 ' 



(AlO) 



(All) 



(A12) 



which is a sort of shifted thermal state and is identi- 
cal to the state obtained in the case of the parametric 
amplifier m% with a mean number of photons given by 
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FIG. 1. Scheme of the experimental apparatus required for the generation and detection of entangled superpositions of 
macroscopically distinguishable states: The idler beam (fc2) of the first nonlinear crystal NLl is used to inject a second 
nonlinear crystal NL2, while its signal beam (fei) triggers the photodetector Di after passing through a polarizer. Modes k2 
and fca are placed within couples of mirrors. The detection apparatus — a rotator, a polarizing beam splitter PBS and the two 
detectors Dc and Dd — probes the field along k2, partially leaking through one of the mirrors of the parametric oscillator. 
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FIG. 2. Photon-number probability distributions for a pho- 
todetection experiment on mode 2e. Pii{n2e), Pv{n2e), and 
P{n2e) [see Eq. (3.5)] are plotted, respectively, in (a), (b), 
and (c). 
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FIG. 3. Scheme of the experimental arrangement needed 
for the measurement of the Wigner function (see text). ^2 
and ks represent the cavity modes of Fig. |j 
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FIG. 5. Tridimensional plot of the Wigner function for 
the d+ mode for an initial mean number of photons 
(n) = N = 2.93. In (a) and (b) two different viewpoints 
have been chosen to display more clearly the quantum super- 
position character of the Schrodinger-cat state. 
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FIG. 4. Tridimensional plot of the Wigner function for 
the d+ mode for an initial mean number of photons 
{n) = N = 14.94. Two different viewpoints have been chosen 
to display more clearly the quantum superposition character 
of the Schrodinger-cat state. 
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FIG. 6. Probability distributions P{x) for the quadrature 
operator Xd_^ ~ (d+ + d+)/\/2 of the d+ mode and for three 
values oi (n) — N. 
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FIG. 7. Probability distributions P{y) for the quadrature 
operator yd^ = i(d^ — d+)/y/2 of the d+ mode and for three 
values oi (n) = N. 
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